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I Abstract 

. We consider the asymptotic evolution of a relativistic spin-^ particle, i.e. a particle 

whose wavefunction satisfies the Dirac equation with external static potential. We prove 
that the probability for the particle crossing a (detector) surface converges to the probabil- 
. ity, that the direction of the momentum of the particle lies within the solid angle defined 

by the (detector) surface, as the distance of the surface goes to infinity. This generalizes 
, earlier non relativistic results, known as flux across surfaces theorems, to the relativistic 

■ regime. 

o 

O ■ 1 Introduction 



In scattering experiments the scattered particles are measured at a macroscopic distance, but 
the computations of scattering cross sections are based on the distribution of the wavefunction 
in momentum space. Therefore a relationship between the crossing probability through a far 
' distant detector surface and the shape of the wavefunction in momentum space is needed. 

This relationship is given by the flux-across-surfaces theorem, which - as a problem in 
mathematical physics - has been formulated by Combes, Newton and Shtokhamer see also 
p. For scattering states (material on scattering states for the Dirac equation is in [Q) 
the theorem asserts that the probability of crossing a far distant surface (physical interaction 
with the detector is neglected) subtended by a solid angle is equal to the probability that the 
scattered particle will, in the distant future, have a momentum, whose direction lies in that 
same solid angle. Moreover, the probability, that the particle will cross the detector within a 
certain area is given by the integral of the flux over that area and time. This has been proven 
for Schrodinger evolutions in great generality, see for instance [|[ ^ 0, ^, ^, 

We consider here wavefunctions V't G L^(M'^) which satisfy the Dirac-equation (conve- 

niently setting c = h = 1) 

^ -i^aidiil;t+ 4-i^t+ Pm%l>t = H-ipt (1) 



dt 



where 



oj'^=lo _J;^ = 1'2'3 (2) 



ai being the Pauli matrices: 



,10 1 I 1 \ -i 

-1/ Vl 0/ \i 
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1 the 2 X 2-unit matrix and 4 the 4-potential in the form 

4 := ^0 + A • a 

with OL := (ai, a2, a^) . In the following we will always denote solutions of the Dirac equation 
by "0* and by ipo the "time zero" wavefunction. 

4 is an external static four-potential, which satisfies condition A (see which concerns 
smoothness and is for the sake of simplicity taken stronger than needed: 

Condition A 4(x) € C°° 3M,C>0: | 4(x) |< M(a;)4+« . (3) 



The norm | • | is defined as: 

\B\:= sup \\B^\\, 

where 

II ^ IU:= 

with the inner product in spin space 

4 

(., •) : C II ^ II. (x)4 (g) ^ C x) E • 

1=1 

Often we have spinors depending on x, in that case we have || V' lis (x). 

The continuity equation involving the quantum flux of a relativistic spin i particle reads 

^V^.V't^V-j, (4) 
whereas the 4- flux is defined for any (p E i^(R'^) C' by 

J= (^") , (5) 



with a = 

For notational convenience we sometimes omit the dependence on x. Furthermore we have 
the usual L^-Norm on the space of 4— spinors given by 



\\v\\=[J W^Wid'x 

We introduce the Fouriertransform of ip{'x.) as representation in the generalized basis ( p^ of 
the free Hamiltonian ,i.e. 

^,(k)= /(27r)-i(^^(x),^(x))d3x ^(k) :=E4<^,(k). (6) 

s=l 

We denote by x the euclidian length of x. 
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We assume that asymptotic completeness holds, i.e. that the wave operators exist on the 
spectral subspace Tiac of the continuous positive spectrum ("scattering state") of the Dirac- 
Hamiltonian: Let i/'out denote the wavefunction of the free asymptotic of a scattering state ■0 
then 

hm II e-*^«Vout - e'^*0 ||= . 

t — ^CJO 

V'out is given by the wave operator: 

Q.+ = lim e'"*e-'"°* ip = f^+V'out ■ 

t — >oo 

The existence of the wave operators and asymptotic completeness has been proven for short 
range potentials. See e.g. Thaller |Q. 

We remark (see Lemma |3.4| (d))(^o[), that the Fourier transform ipout,sO<) of V'out equals the 
generalized Fourier transform tpl of i(j in the generalized eigen-basis of the Dirac hamiltonian 
with potential. 

In general, we do not have much information about scattering states. One can prove the flux 
across surface theorem with conditions merely on the "out" -states, where the corresponding 
properties of the scattering states are hidden in the mapping properties of the wave operators, 
or, better, in the smoothness properties of the generalized eigenfunctions. On the other hand, 
one would like to be sure, that such conditions are not too restrictive on the set of scattering 
states. 

We introduce the set Q of functions ipout, for which the flux across surfaces can naturally be 
proven: 



f(^<^^G^^ 3MeR: || a^/(k) ||,< Af(fc)-" forj = 0,1,2; neN 

•'^^ ^ \ VkT^O: II ^'^'-^^^/(k) ||,< Af(fc)-" forneN ^ ' 

where 7 = (71; 72; 73) is a multi-index with| 7 |< 2, :— d']}^d1^d1\ 

This set maps under the wave operator to a dense set in the set of scattering states. After 
the theorem we shall give under more restrictive conditions more detailed information on the 
set of scattering states for which the theorem holds. 

The paper is organized as follows: In the next section we shall state the theorem. We shall 
also give its formulation in covariant form, but we shall prove the theorem using the rest frame 
of the dectector and the potential. 

The following sections contain the proof of the theorem. We first prove the statement for 
the free case (4 — 0) and then for the case of nonzero potential. Both are done in section 3. 
The proof relies almost entirely on the stat iona ry phase method, which we need to adapt to 



Appendix 4.1 



our purpo ses. The main lemma is lemma (3^), whose lengthy technical proof is put in the 



The difficulty we have to face and which makes this paper not a simple generalization of the 
results in the Schrodinger situation is, that the time evolution with the Dirac hamiltonian is 
not of a "nice" form for the stationary phase method to be easily applied to. The Schrodinger 
case is easier. On the other hand, the expression for the flux needs no differentiability of the 
wavefunction and one might be lead to believe, that to describe scattering in the relativistic 
regime is simpler — in particular less restrictive theorems should result. One may even get the 
idea, that asymptotic completeness and the flux across surfaces theorem become more or less 
equivalent statements in the relativistic regime. But we are far from that. Nevertheless, that 
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we require smoothness and good decay on the potential may well be due to our method of 
proof. 

We also need information about the generalized eigenfunctions of the Dirac hamiltonian 



with external potential, see Lemma 3.4, whose proof is also put into the Appendix 4.3. The 



appendix, which in fact is almost half of the paper, contains other tedious technical details. 

Acknowledgement: Very helpful discussions with Stefan Teufel are gratefully acknowledged. 
We would also like to thank the referee for a very critical reading of the manuscript, which let 
to a clear improvement of the paper. 



2 The Theorem 

The flux-across-surfaces theorem deals with the flux j integrated over a spherical surface at a 
far distance and asserts that 

1. the absolute value of the flux and the flux itself yield the same asymptotics, allowing to 
interpret the flux integral as crossing probability |, |ll|, 

2. the crossing probability equals the probability for the momentum to lie within the cone 
defined by the surface. 

Theorem 2.1 Let ^ be a scattering state with outgoing free asymptotic tpout, whose Fourier 
transform Tpout lies in Q(cf. Let R^dil be the surface element at distance R with solid angle 
differential dQ and let n denote the outward normal of the surface element. Furthermore let S 
be a subset of the unit sphere. Then for all ti £ R; 

lim / / j(R,t)dtR^dn = lim / / i{R,t) ■ ndtR^dn 
Js Jti Js Jti 

{iioutO^),iiout{k))Pdkdn , (8) 



s Jo 



Observing that || "tpout || (k) does not depend on time, we can choose a coordinate system 
t' = t — ti, so that we may for definiteness always put ti — in (^) . 

The conditions on "fAout can be translated into more detailed conditions on the scattering 
states under more restrictive conditions on the potential: Let 



Condition B | 9"4(x) |e L^{R^) Vn e {0,1,2...} 3M|4(x)| < Mi 



-6 



Then (for the proof see Appendix 4.4) 
Lemma 2.2 

$out{k) eg^ V(x) e g (9) 

where g is the space of functions ipi^) G T^ac with x^y '""ijj{'x) G L^ for all j = 0, 1, 2; n G No, 
where y :— ~iJ2i=i'^i^i- 



4 



2.1 Covariant form of the theorem 

As we deal with a relativistic regime, it might be of interest to have also a covariant formulation 
of the theorem. As (V'oiit, V'out) is not conserved under Lorentz function we use 

of which it is known, that {^out^^out) is a Lorentz-scalar (see for instance Then the 

flux-across-surfaces theorem reads in a general and covariant way: 



Theorem 2.3 Let the conditions of Theorem 2.1 he satisfied. Let 
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xoy:= xoyo -'^Xjyj 



be the Minkowski scalar product. Then for any subspace Z C {x \ xo x = m } C M"* and any 
smooth scalar function r/{x) bounded away from zero: 



lim / j(x)onda= / (^i^,(fc),V'i',(fc))da. (10) 



where 

Z{X) := {y\3x€ Z ■.y = Xt]{x)x} C 

and da is the invariant measure on Z, da the invariant measure on Z and n is the vector 
orthogonal on Z with Lorentz length one. 

This formulation may perhaps not be directly guessed, but once one understands its basics 
like (p^, this formulation becomes clear: The arbitrariness of the scalar function rj follows 
directly from (^8|), observing that 

lim i^{Xk) = lim ip{Xri{k)k). 

X — *oo A — *oo 

Physically this is related to the fact, that (on big scales) it is possible to "catch" any part 
of the wave-function in different ways (for example by using a detector which is "close" and 
catches the wavefunction at an " early" time or one uses a far detector at a later time-interval) . 

Let us explain how (H) follows from (|l0|). We choose a set Z whose projection on the 
t = 0-subspace is a cone with angular distribution S: 

Z = {k\ ^ € S}r\{kok = mH . 
k 

The invariant measure on the mass hyperboloid da — ^^^^^^ we get for the right hand side 
of (0) 

' {i^il,{k),i,^l,{k))d^ ^ I I {$out{k)Jjontm k^dkdQ . (11) 



L 



z Js Jo 



For the left hand side of (10) we take: 



r]{x) :— — X ^ 

X 
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As both integrands in ( pX)| ) are bounded, a small neighborhood of x = can be neglected. For 
constant A, Z represents a radial surface with arbitrary time i > 0. So we have: 

lim / j(x)ond^= lim / / j(Il,t)dtR'^dn . (12) 

3 The proof 

3.1 Scattering into cones heuristics 

The flux-across-surfaces theorem is based on an asymptotic connection between the shape of 
the wavefunction in momentum space and in ordinary space. This is often referred to as the 
scattering into cones theorem, which has been proven for non-relativistic particles by Dollard 
| p2[ . For that one chooses a certain parameterization of and evaluates the wavefunction, as 
the parameter of the parameterizations goes to infinity. In the non-relativistic case, it is easiest 
to choose time as the parameter of the parameterization. In the relativistic case it is simplest 
to have lorentz-invariant three-dimensional subspaces of the time-like part of as leaves of 
the parametrizatiori|^ This can easily be done, by choosing a lorentz-vector as argument of 
■0, i.e. a vector x with a;oa; = a::Q — x-x = Am^. Set ipi^k) = = Ak, t = \\/k? + m?). 
We denote the two different eigenstates of momentum k of the free Hamiltonian with positive 
energy by (^^, whereas the s labels the two different spins our electron may have. In the 
standard representation these eigenstates can be written as: 

(Pk = e Sk. (13) 

where the ^'i'^- 



4 - {2EkEky^ 





fcl 



4 = {2EkEk)- 



where 



k — k-} ± ik-i 



El. = El, + m 



Ek = ^/B + to2 . 



(For a detailed calculation of these spinors see (jj])) 

The asymptotics result from a stationary phase analysis: 



V'(Afc) ^ J7(t = AVfc2+^)V'(Ak, 0) 



Y^^-rHx^k-+m- I (27r)-t(p^,(Ak)V;,(k')d3fc' 

s=l 
2 



(27r)-Se''''-^''4,V^,(k')d^A:' 



For convenience we define: 



^We only parameterize the time-like region, as for big time-scales the main part of our wavefunction will be 
in this region. 
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This leads to: 
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In view of the stationary phase method, in the hmit X oo only a small neighborhood of the 
stationary point of the phase function 



/i(k') := (VA:'2 + mVfc^ + m2-k' • k) 
will be relevant for the integral. The stationary point is given by: 

Vk'ft.(kstat) = =^ kstat = k (15) 

Without loss of generality we can set k2 — — 0. Near the stationary point the phase is to 
second order: 



-iX{Vk'^ + m^Vk^+rn^ - k' • k) « -^(m^ + J"' (fc^ - kf + \{k'i + k'^)) 



This in equation (|l 41) leads to 

^{Xk) « j (2^)-ie"'^'™'+^(^^^'^"'^' + ^^'"+'"^^V^(k')d3fc' , 
and replacing "(/"(k') by "(/"(k) we obtain by integrating the gaussian 



-i\m' _ j 



V'(Afc)«^-— ^^(k)W^ + l 



(zA)3 



m 



We shall state now the stationary phase result in a some what more general setting, to cover 
also applications to the potential case considered later: 

3.2 The stationary phase 

Lemma 3.1 Let x be in Q (see^) and let the "phase function" g be 



g{k') = y/k/^ + m2 + a | fc' | -y • k'. 
Let kstat be the stationary point of the phase- function: 

Vg{kstat)=0. 
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Then there exist Ci, C2, C3 G M so, that for all x with \\ dlx \\s<\\ d{x \\s for j = 0, 1, 2, y G 
and a > 

II / e-^'^s(''')x(k')d'fc' - CiM-^x(k.tat) ||.< C^lJi-^ + C^i^xiKtat) ■ (16) 

_For phase functions without stationary point Ci — C3, — Q. Moreover the Cj are uniformly 

bounded for all x, a cind y. For a ~ we can choose Ci = (— 27ri)5e~*^^(''''*"')-^^^^^^!^-^!^^^^-^ and 
C3 = . 

One may be disturbed about the nature of the inequahty ( |l6|) when C3 ^ 0. The point here 
is that our estimate is uniform in fc,/„/ and then later we shall use for C 7^ such that 
kgtat will be of order so for C3 7^ the last term will be part of the leading term in our 
estimation. 



This statement is a sli ght adaptation to our situation of a theorem of Hormander 1 14 , and 
its proof in the appendix 



4.1 



3.3 Scattering into cones for a free particle 



Applying Lemma 3.1 to (n4j) we choose: 



= XVk^+m^- a - 0; y - ^ ; x(k') = (2^)-^^(k') 

and calculate the stationary point kstat- 

—^^^^=-y = 

kstat = y^ik%at+m^) 
ym 



obtaining 

Corollary 3.2 ("Scattering into cones") There exists a constant C < 00 so that for all k G M'^ 



II V-CAfc) - + 1 IU< 



Note, that this implies 



3 ^ 

lim supdl %/A ^(Xk) \\s - II ^ik)J— + 1 lis) - . (17) 
A— ►oo k V 

For the flux-across-surfaces theorem we need the asymptotics of the relativistic quantum 
flux (ph of the particle. Since all the ai are bounded matrices and ijj G G, we obtain from ^ 
and (O) for the flux: 
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lim sup I X^jiiXk) - {^{■k),aii^ik))i^ + 1) 1= . (18) 



Next observe (see the appendix 4.2), that: 



(^^(k),aV^(k)) - -=L=(v7(k),v7(k)) . (19) 
V fc + 

Thus we get the uniform bound: 
Corollary 3.3 

Ve > 3A e M : 

sup I X^UXk) - (^(k),V^(k)>^Vfc2 +m2 |< e . (20) 
k 

Observe, that after a long time of propagation, the flux at x = Ak will always be parallel to k. 
So in the limit i — > oo it will always point away from the origin of the coordinate system. 

3.4 Flux across surfaces for a free particle 



Theorem 2.1 reads in this case 



lim / / j{Ti,t) -ndtR^dn- / (V'(k), i/'(k))fc^rffcdri = (21) 
Js Jo Js Jo 



and 



R — ^oo 



SJO JsJo 



lim / / jCR,t)dtR^dn- / {'^P(]^),^P{^k))k^dkdn^O . (22) 



In the following, we will prove (|2^) by inserting the longtime asymptotic ( PQ ) for j and showing, 
that the integral of the error we get by this approximation tends to zero in the limit R oo. 

Now, the long time asymptotic of j is parallel to the normal n of the radial surface. Therefore 
the longtime asymptotic of j is equal to the longtime asymptotic of j • n. More detailed, one 



sees that using the approximation ( pO| ) for j in ( |2l| ) and (22), the bound on the error terms in 
( pl| ) and (|2^) arising from ( pO| ) are equal. 

So the proof of (^2|) is essentially the same as for (^l]) and we shall concentrate only on 
showing (pi]). 

The left side of ( plj ) includes an integral over t, whereas the right hand side is integrated 
over k. We therefore substitute for t in the flrst term, to get integration over k, too. Since A 
plays the role of a time parameter it is natural to substitute: 

Rn 

with 

A — . 
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But this substitution is only possible in the time-like region (t > R). So we first handle 
the integral starting at t = R, later we deal with the space-like part of the integral. Then, 
substituting t by k, we obtain 



S JR 



i{R,t)-ndtR'^dn = / / i{R,-Vk^+m^)-n 
JsJo k 



-dkdn 



s Jo 



2 

{\{k)k, \{k)y/k^+m^) ■ n kX{k)^dkdn 



The integrand is now in the form that we can replace it by the asymptotic in (|20|). 

It turns out however, that the error in the integrand will be ~ ,,J', ^ which is not inte- 
grable, therefore the replacement is not straight forward. We separate large momenta k > X 
and small momenta k < X. In the following we choose X > m. Given X and Rq — XqX 

k<X^^=X{k)>Xo = ^. 

k X 



Then by (gO|) for smaU momenta {k < X ^ t > R^ 1 + 

Ve > 3i?o e R Vi? > Ro 



'j(R,t) -ndtR^ - (^^(k),V(k))fc2)dfcdr2 I (23) 
Jo 



rX 2 ^ ^ 

j(Afc, XVk^ + m^) ■ - (V'(k), ilj(k))k^dkdn 



where 



Is Jo VP +771" 

< [ [ /""'^ dkdn^:x{X)e 
Js Jo V fc^ + 



Jq VK + TTT- 



Given X we can take e arbitrarily small, choosing Ro large enough, so that the r.h.s. of (23) 
goes to zero. Thus 



lim lim I / / j{R,t) -ndtR^dn- / {il^(k),i^{k))k'^dkdil \= . (24) 

' S J rJ 1+^ JSJO 



X — >oo R — >oo 



For the large momenta note that by virtue oi tj} ^ Q: 



X 

and all it remains to show is that 



lim / / (7A(k),7A(k))fc2dfcdl7 (25) 
Js Jx 



lim lim / / j(R,i) • ndti?^dl7 (26) 



X — *oo R — *C30 



s Jo 
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where we also included the time integration outside the Hght cone, which we excluded in the 
substitution. 

We first estimate the part of the integral (|2^) that lies in the space-like region (more precisely: 



t e [0, R]) then we estimate the time-like part near the light cone ( t G [i?, Ry 1 + yr])- That 
is, we first show that 

lim / / i{R,t) ■ndtR'^dn = . (27) 

That this holds is physically related to the fact, that a particle moves slower than light, so 
for big time and space scales the main part of the wavefunction will be inside the light cone. 
This follows from a straightforward application of the stationary phase method, outside of the 
stationary point. Two partial integrations lead to: 



i^(x,77x)|U = II /(2^)-ie-"(^'^^''-^-^)V^(k)d3fc|j, 
(27r)-ie-*^f^?(k)d3/fc ||, 



where 



Since 



it follows: 



< — / 27r) 2 — + — )\^d-'k 

a; J g 5 9 9 



, kiT] I fci I 

-5=1- > 1 - , > 



II V(x,r;x)|U<(27r)-i^ (28) 



uniform in 77 < 1. Hence 



R 

Jim / / i{Ii,t) -ndtR^dn 



< 



4n lim / II ^(x, t) II 2 dtR' < -^C| lim R'^=0 
It is left to prove that the second part of the integral in (Eq) goes to zero, i.e. that 



X — *oo R — *oo 

The scalar norm of t) is: 



lim lim / / i{R,t) ■ ndtR^dfl = 

s JR. 
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H(x,t)||, = II / (27r)-ie-'^^^*+*-^^;d3fc lU (29) 



II / (27r)-ie-'(^/'='+'"'-k.r)i^^3^ |j^ (30) 



Applying Lemma 3.1 with 

^l = t■ a = 0; y = v- x(k') = (2^)-i^(k') 

we have by (p^), that: 

g-»ij.t+^k.x^(^)^3;. - |U< Cat 

As '(/' is bounded, we have: 

3M eM:yt> R II i^(x,t) ||, = || y e-*^^*+*-"V^(k)rf 



So 



So we can write: 



/ j(R,t) • nR^dn \< An 
Js 



MR^ 



/ / 3{R,t) -ndtR^dn 

Js Jr 



m 



2 



< 2tiMR\R-' - R-'[l + ^y') = 27rAf(l - (1 + ^)"') • 
This term goes to zero as X ^ oo 

3.5 The flux-across-surfaces theorem with potential 

3.5.1 Generalized Eigenfunctions for the Dirac equation with potential 

For the proof of the free flux-across-surfaces theorem we used the tp^ as basis of the Hilbert 
space. In the potential case we adopt a new basis for doing calculations. 

Like in the free case, we again get four linear independent eigenfunctions for each k, two 
of them have positive energy-eigenvalue E"^^ = Ek — Vfc^ + w^, two of them have negative 
energy-eigenvalue E^^^ — —Ek- We denote by <^^(x) eigenfunctions with s G {1, 2}: 

EkKi^) ^ {Ho + 4M{^) ■ (31) 
The corresponding Lipmann Schwinger equation reads: 

^^(x) = ^^(x) + {Ek - Ho)-'4K{^) ■ (32) 
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We replace the formal expression {Ek — Hq) ^ by the integral kernel G^: 

(ii;^ - i/o)G+ (x - x') = <5(x - x') . (33) 
The expUcit form for G^(x — x') can be found in 



3 3 



ikx 



G+(x) = i-e^'=- +^«^.fc^ +^^) ^: ^5+(x) . (34) 

Thus: 

^^(x) = ^^(x) - I 4(x')G+(x - x')^^(x')d3x' . (35) 
For S^, defined in (^), we have: 

3 3 



for j = 0, 1, 2. Choosing a; > 1 we have 



— < 1 ^ < 1 



and it follows, that 



Thus with 



we have: 



dlS+ \<\ + ^ a,{k + 1) + Pm) I 

1 ^ 



47r 



I 9^:5+ |<| 9^5+ I (37) 

for j = 0,1,2, x> 1. 

For the next steps we need some properties of the generalized eigenfunctions . Wc summarize 



these properties in the following Lemma which is proven in the Appendix 4.3 



Lemma 3.4 Let ^ satisfy Condition A (Qj. Then there exist unique solutions ^l^i^) of ^5 
for all k gM.^ , such that: 

(a) For any k G M"^, s = 1,2 the functions ^^{'x.) are Holder continuous of degree 1 in x. 

(b) Any ^^(x) which is a solution of (pqj automatically satisfies i\3^). 
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(c) The functions 

a(x) := ^^(x) - ^£(x) (38) 

are infinitely often continuously differentiable with respect to k, furthermore we have for 
j £ N and any multi-index 7 with | 7 | < 2 ; 



i) sup II a;Ck(x) ||s< 00 
^i) sup II 9^ ^l^f |U<oo 

XGM3 I X + 1 |J 1 



) sup II fcl7|-lj^7 aw ^ _ 



(d) T/ie form a basis of the space of scattering states, i.e. for scattering states '4>{x,t): 



2 . 

V'(x,t) =Y,J {27T)-h-'^^^'^U^)$outA^)d'k (39) 



V'o«t..(k)= y (27r)-5(^^(x)Xx))d3a; (40) 
where V'ojit,s(k) is t/ie fourier transform ofipout — f^+V'- 
3.5.2 Flux-across-surfaces for the Dirac-equation with potential 



We prove now Theorem As in the free case only the equality of the second and third 
integral is shown. From the nature of the estimates in the proof it will become evident, that 
essentially by the same argument as in the free case, the first equality can be established, and 
we do not say anything more to that. 

We again split our flux integral into two parts, one inside the light-cone (from R to 00) and 
one outside the light-cone (from to i?), where the main contribution comes from the times 
t > R, i.e. we prove that 



00 



i) lim I / / j(Ti,t)-ndtR'dn- / (V'out(k), ^out(k))fc^dfcdri h 



SJR JSJO 
R 

2, 



ii) lim / / j{Ii,t) ■ndtR^dn = (41) 



We start with i): 
According to (p9|) 



s=l 

Setting 



2 

^(x,t) = J2J (2^)-ie-'^^^*^^(x)^„,,,,(k)d3fc . 



V'out(k') = ^4,?Ao„f,5(k') 
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and using (|35|) with ( |38| ) we get: 

=: 5o- 51-^2. (42) 

5o is the propagation of the free outgoing state. The free Flux-Across-Surfaces-Theorem yields 
therefore: 

hm I / / {So,aSo)-ndtR^dn- f f (^out(k)>out(k))fc2dfcdri h . 



R — >oo 



SJR JSJO 



Hence for (i) it remains to show, that (using ||) : 

lim / / {i{R,t) - (So^aSo)) -ndtR^dn 
Js Jr 



■s-'fi j=0 j=0 



= lim / / {{iP,ay S,) + {y S,,ail;))-ivdtR^dn = Q . 
'^^^JsJr f-{ ^ 

By Schwartz-inequality we need only show: 

r 

lim / / II V lU y II 5*^ II, dtR^dn ^ . (43) 
R^oo Js Jr ~^ 

We first want to estimate || 5*1 ||s. Recalling we get by Fubinis theorem: 

Si = y(27r)-ie-^^^^*y^J^^^5+(x-x')4(x')e"""'d'x'^out(k)d3fc 

(2^)-ie-'^fe'+'"'* ""'^'^ " ' g+(x - x')4(x')e*'^"''d3xVout(k)d^A: 
I X — X I 

=: /(27r)-i ^ gi(x,x04(xOrfV 
J I X - x' I 

where 

5i(x,x') = J (27r)-ie-*^'='+"'*e*'=l''-^'IS'+(x - yi')e*^-''' d^x'$out{k)d^k . (44) 
Next we use Lemma ^.ll , setting: 
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^ = i; a = X - x' I; y = t^^x'; k' = k; x(k') = (27r)-5S'+(x - x')V'(k') . 
With regard to , the function 

x(k) = (2^)-i5+^(k') 

satisfies the properties we need in (p^). Furthermore we observe that for the stationary point: 

h . . 

a — y = 



kstat = ^Jk'^^t^t +m'^{y - a) 

So we can estimate kstat by: 



fcLt + "^'i"'(^'- |x-x' I) < ^kl,, + m^xt-' . (45) 

Hence by we obtain for ( ^4|) that there exists Afi < oo, bounding in particular 

\/kstg^t + ™^x(kstat)i which is bounded by the choice of V'out G 5 and incorporating also the 
constants Ci and C2, uniformly in y and a so that: 

II ^1 II. < II Aht-^i + x^) f - — L^4{x')d'x' lu 

J I X - x' I 

= Afii-iPi(x) ^™ . (46) 

That the function Pi goes to zero in the limit a; — > 00 may be seen as follows: 
For any function /(x) G with limsup^^^^ | x^/(x) |< 00 we have: 

hm X I / ^-—f{^')d^x' I 

x^oo J I X — x' I 

< Mm X f \ ^/(x - x') I d^x' 

x^oo J X 

= lim x{f I l/(x-x') |dV+ / I l/(x-x') |d=^x') 

< hm :r(sup{|/(x)|} / l.d^^' + -f | /(x - x') | d^x') 

< lim ix^ sup{| fix) 1} + lim 2 / | /(x - x') | d^x' < 00 (47) 



where i?(a, r) means the ball with center a and radius r. 

Next we estimate || 5*2 ||s- According to ( ^ we can write it down as: 

^2 = E / (27r)-ie-'^^^* J ^^^^5+(x - x')(.x' + l)4(x')|^d^xVo„*,,(k)d3fc 
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Therefore we again use Lemma 3.1, setting: 



^ = t; a = X - x' I); y^O-k'^k; x(k') = (2^)"i ^ f^^k " x')^^o„M(k') 



With regard to ( |37|) and Lemma 3^(c) there exists a M2 < 00, so that the function 



X=(2^)-^M2^+^(k') 

satisfies the properties we need in (pj|). 

Since our phase function has no stationary point we get with (|l|): 



^2 \\s<M2t-^ 



J I X - x' I 



Choosing (x' + l)4(x') for / in the most left side of ([47|), one can see, that a;P2(x) is bounded, so 
P2 goes to zero in the limit a; — > 00. Since So is the analogue of the freely evolving wavefunction, 
we have by Corollary |3.2|: 



We use the estimates 
M0+M1 + M2: 



\\So\\s<Mot-^ . (48) 
and (Esl) in the right side of (E^ and get, defining M :— 



lim I / / lUII \\,)dtR^dn I 

R^°^' Js Jr. ]-[ 

< lim / M2(Pi(R) + P2{Ii))t-^dtR^ < lim 3M^{Pi{R) + PalR)) = 

R — >oo J R — >oo 

and ( ^ ) is proved. 

Like in the free case, (pd| ii) follows directly from an analogous argument which used equation 



8|), thus we prove (ES) for the case at hand. Since in (41 ii) we need only estimates of the 



wavefunction for times t < x we have in view of ([l^, setting t — rjx with < 77 < 1 and using 
Fubinis Theorem: 



V'(x,r?a;) 



(27r)-5e-*^'='+'"'''^'+*'^-''^out(k)d3fc 

-zVfc^+m^?7a+ifc|x-x'|+ik-x' ^CxQ-Si^lx - x')^out(k) 

(27r)i I x-x' I 
(x05+(x-x 
(27r)i I x-x' I 



d^kd^x 

d^kd^x 



s=l 

= : S'o — 5*1 — 5*2 . 
For Sq we have (|2^), for the other summands we define 



-iVfc^+rra2»7a: + tfc|x-x'| 4(xOCk (xQ ■S'k (x - x')V'o„^g(k) ,3 ,^^3^, 



^1 := /(27r)-5e 



5+(x-x')V^out(k)d3A: 



2 

52 := ^ / (2^)-ie-^^'='+™'''"+*'=l"""'l+''*""'e-*-"'Ck(x')^k (x - x')V^o«t,.(k)d3A: 

s=l 
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So we have for Sj, j=l;2: 



s, 



We can estimate the Sj by two partial integrations. One can easily see, that the phase functions 
of Sj have no stationary point. This leads to: 



where 



Since 



(27r) 



^'('^)x,(x,x',k)d3fc 11, 
II /(27r)-Se-"s('')9fc,(iafe,^)d3fc||, 



9 9 



ff(k) 
Xi(x,x',k) 

X2(x,x',k) 



X — X 



5+(x-x')V^out(k) 



2 

s=l 



x' fcl I X — x' 

a:: fcx 



2 - fc ^ a; 



fc 



Vfc^ + TO- 



=)>o 



is bounded and due to Lemma 3.4 the Xj bounded, we can find C2 < 00 with: 
So is bounded (see ^) and the analogue of ( p8| ) is proved. 
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4 Appendix 



4.1 Proof of Lemma( |3.1| ) 

We consider for a family of phase functions g, which we should think of being indexed by 
a> 0,y: 



f;(k) = \/k'^ + rn? + a | fc | — y • k 

the integral 

/ := y e-*^9(k)^(j^)^3^ 

where x & G (seeQ). 

We shall find its asymptotic behavior as a function of fi. In major parts we will recall 
the proof of theorem 7.7.5 in the book of Hormander ]T^ , which unfortunately is formulated 
for compactly supported x ^nd which moreover does not give uniformity over the family, 
i.e. uniformity in a,y which we need. The compactness can easily be handled but for the 
uniformity we must invoke the special form of the family of phase functions g and we shall give 
the argument here. 

The stationary points of the phase functions are given by: 



//I \ ^stat , ^stat ^ 

gi^stat) = n +^T y = 

VKtat + "1 

_ m{y-a) 

Ktat II y. (49) 

Since kstat is a function of a and y, we sometimes use the phrase: uniform in k^tat to express 
uniformity in a and y. (I) For y > a+l there is no stationary point and for y = a the stationary 



point is at k^tat = 0. 

First we handle the family where y £ [a + ^; a + 1[. These phase- functions do exactly have 
one stationary point bounded away from zero: 

Later we will handle phase functions, where the stationary point is close to zero and phase- 
functions without stationary point. 

We choose a coordinate system, where the fci-direction is parallel to y. So the stationary 
points will have the coordinates {kstat, 0,0). To estimate the integral, we separate from the 
integral the contribution coming from near the stationary point. This part of integral includes 
the leading term. Therefore we define a smooth function pk^tat which is one near the stationary 
points and zero away from the stationary point. (We shall omit further on for ease of notation 
the index kstat)- 

More precisely we define the compact set Q by: 
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keQ^fci e [^,2k,tat]Ak2,k3 e [-1,1] 

and choose 

p(k) := 1 Vk e Q (51) 

falling quickly off to zero outside of Q, lets say 

/)(k) := VA; ^ (52) 

where is some e-neigliborhood of Q for some e > 0. With the help of p we can split 
X = Xi + X2 by defining: 

Xi '~ p\ X2 (1 — p)x 

7i := Je-'^'s(•^'>XlO^)d^k h := / e-^''s(k)^2(k)rf3fc . 

This split has the following advantages: 

The compactly supported xi includes the stationary point, so /i can be estimated the same 
way as in Hormanders theorem, but with focus on the uniformity of the estimates. X2 is zero 
near the stationary point, so I2 can be easily estimated by partial integrations, p has been 
defined in such a way, that we may estimate the terms we get by the partial integrations 
uniform in kstat- 

We start with Ii. We move the stationary point to the center of our coordinate system 
setting k' := k — k^tat, i-e. (/(k) becomes 5(k') = g{\i! + \istat)- Slightly abusing notation we 
simply write ^(k') for g. By Taylor's formula we obtain a function /: 

5(k') = 5(k' = 0) + ^ "^"^'^^^ = Q^'^"' + m , (53) 

l7l=2 

where bounded. 

Computing the second-order terms of 5((k') we find that only diagonal terms survive at 
(/Jstat, 0,0) and 



so that 



We define: 



a^,5(k' = 0) = dlg{\^^\^stat) 



fc^ - fc? + TV? fc2 - A:? 

(— ==^ + W=^s... (54) 



dl,g{\^ = Q) = ^'""+"'' 3 +ar^forj^2,3 

a,V(k' = 0) = — =^^. (55) 



E;Lia,2'5(k' = 0)fcf 
52(^, e) := ^ ^^^^ . (56) 
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By this definition, 52 does only depend on the angular, not on the radial coordinate of k'. 



Using ( p6[ ) in (53), we may write 

1 



g(k')-g(O) + ^fc%(.9,0)+/(k') . (57) 
Furthermore for s G [0, 1] set: 

gs :=.g(O) + ifc'2g2(^,0) + s/(k') (58) 

and 

Note that g = gi, Ii = /(I). By Taylor's Formula there exits ^ < 1 so that: 

/i =/(l) = /(0)+5,/(s) 1^ . (59) 

We begin with /(O), introducing spherical coordinates. With slight abuse of notation: (leaving 
the notation for the functions unchanged) 

/(O) =y"e-*'^(^'(")+^'="f-(^'«))xi(fc',z9,0)fc'2dfc'dr! . 
Writing xi = x{k' = 0) + x the integral splits into: 



/(o) = J e-^'^(f(")+^'="5^(^'«))x('fc' = 0)k'^dk'dn 

+ J e-^M(s(o)+ifc'^s.(tf,e))~(^/^ e)k'^dk'dn =: + if . (60) 
The integral is a gaussian integral, which includes the leading term: 



II = / e-^''(»(°)+^'= ^^(''^^"^(fc' = 0)k'^dk'dn 

-^^llU ^^.f «('^'=o))'=?,-...(o)^(k' = 0)k'^d^k' 



(2^)^M-^e-''^(°)(na,2,5(k' = 0))-^x(ks*at)- (61) 



For a = the 9j./2(7(k' = 0) terms can be easily calculated. We get: 
dlg{\^ = Q) = dl,g{k^kstat) 

J J 

_k^ + m^- k] 

So we get: 
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ll is the leading term of our integral. For a — we get the desired value for Ci (3.1) 
For /J put: 

which is bounded and smooth. 



/2 = J e-^f 9)k'^dk'dn . (62) 

One partial integration leads to: 



-^,hk'\,^i}.e^ dk'Hk', ^, 0)k'^ + 20(fc^ ^, g)fc' ^^^,^o II 
52(1^,6') 



So another partial integration is possible: 



k'{92{i},0)y 



^ ' (52(^,0))^ 

<- , ^fc'>^(fc^^,g)fc- + 2^(fe^^?,g) ,^,,^^ „ 

< II J dk'i {g2{d,0))^ )dkdn\\s . (63) 

With our definition of Q, the support of the integrand increases and g2{'&,0) decreases poly- 
nomially with kstat (see ( |55|) and (|5^)). While the support moves away from the center of our 
coordinate system. But X = X ^ X^stat) and its derivatives decay faster in kstat than any 
power, so we get a constant C uniform in kstat with: 

II < M-'C . 

For Ii it is left to estimate dsl{s) |^: 

dsl{s) J -ifif{k',^,9)e-'''3(('''^'^-'^)xiik',^,0)k'^dk'dn . (64) 
By Taylor's formula we can define: 

f{k', d, 9) := f{k', d, 9)k'-^ g{k', d, 9) := k'~^dk'9i{k' , d, 9) 

and thus: 

dsl{s)\(^= [ -ifif{k',^,9)e-"'3i('''-^-'^)xiik',^,0)k'^dk'dn . (65) 
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On (see below (|52|), g is infinitely often differentiable. So these functions are well defined 
and bounded on Q^. 

To estimate the integral by partial integrations we have to assure, that has only one 
stationary point, which is k^tat = as one easily sees from (|6^). 
By m-. 



9i = .9(fc' = 0) + h'^g2{^, e) + e/(k') = ^.9 + (1 - [gik' = 0) + ifc''g2(i9, 
Looking at 



(66) 



we observe, that 



dig = dl fc'2 - 2k'Ktat cos(i?) + kl^^ + m2 + k'^ - 2k'k,tat cos(i?) + k^^^ - y k') 

Q , fc' - fcsiat C0S(1?) fc' - fcstat C0s(?9) 

Ok' ( — , + a — - ) 



■ - 2fc'fc,t,t cos(7?) + k^^^ + m2 Vfc" - 2fc'fc,t,t cos(i?) + fc^,,, 



(l-COs(z9)2)fc2,„,+m2 , ^ (1-C0s(^)2)fc2^^, 



> 



- 2fc'fc,t,t cos(i9) + + ^k'^ - 2k'k,tat cos{d) + k[ 



2 



And for k e Q^, fci is positive, so the angular component e] — ^, we also have, that 
on also 52 is positive. Since ^ G [0; 1] it follows, that dl,g^ is positive, so d^g^ is strictly 
monotonous on and has only one stationary point. Recalling the definition of 5 (see |65| ) we 
see, that g is bounded away from zero. 

Now we can estimate the integral (p5[). By partial integration 



g{k , u, fc*) 



^ gik\^,0) gik',^,0) ' 



7.,,,.. ^ /(fc^^,g)xl(fc^^,g) ,, , ^ /(fc',^,^)xi(fc',^?,g) 



Setting 



Hence 



This term is similar to (|62|) . The only differences are, that we have tp instead of 4> and 
instead of go- 

So with the same estimate as in (^) we get: 

II dsl(s) \^\\s< fi II J dk (g(k^ d e)Y d^Ws ■ (68) 



dj{s) |5= / e-"'s^'-'''-'^-'''^'i{jik',^,e)k'^dk'dn . 
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This term again has uniform bound in kstat, as its support moves away from the center of the 
coordinate system. So we get a constant C uniform in kstat with: 



II dsHs) k\\s<f^-^C . 

Now we estimate 7(2) (|53|). As this integral includes no stationary point, two partial inte- 
grations are possible without any problem, but we have to assure, that we can estimate the 
factors we get by these partial integrations uniform in kstat- To be able to find an uniform 
estimate, we estimate the areas of x separately. 

So we again split our integral: 

*:ie-B;|fc2|>l ^fcieB;|fe2|<l;|/c3|>l 

= : + + 

where B := [^;2kstat]- 

The integrals and /| we estimate by two partial integrations under the fci-integral. This 
leads to: 



-2 f II Q^2 , „X2.9 „X25 II ,3; 

= M / \\3— + 3^^-3—\\sd'k 

t2 II ^ „-2 f II 1 n X2(k), „ ,3, 



- ||3% + 3^-3^|Ud3fc (69) 

Jki>2k,tat 9 9^ 9-^ 

where g(k) := ^^^^(k); g'(k) := 9^-5 (k) 

At first sight these estimates do not seem to be uniform in a and y. In fact 



<7(k) = -+a-^ 3 



and 

5"(k) 



are bounded on the area of integration. So it is left to show, that we can find functions hj with 
j=l;2, which do not depend on a and y and which is bounded away from zero on 1R'^\Q with 

/ii(k) < .g'(k) /i2(k)<.g(k) 

for all a, y, k. 
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For this we estimate g for ki < ^^j^. As g > 0, it follows, that (see (60)) 



|^(k) \^y- ^ :-a^>^ 



For fci < and by virtue ?/ > a + ^ > 5 • 



For fci > we estimate, using that 11 ~ a . 2^°'"' ~ 0: 

I 5(k) I = y j,^ , ^ - > 2/ - a 



> y-a 



2 + 

kstat , kstat ki 



kstat k\ 



^ fcstat V/C^ + - fci V ^stat + ™^ 



fc'tat(fc'+^')-fc?(fcLt+"^') 



Recalling fc e [0; 



l<?(k)| > 



i'^stat'"' 



V3 



(kstatVk^ + m2 + fci^A:2^„t + to2^ ^A:2 +m2^A:2^„t +771 

3to2 



I 5(k) |> ^""^ , =: /ii . 

8(VFT^ + 2fci) y^fTT^ 



For fci > 2kstat, g' is positive. Therefore similar as before: 



f'WI = /,2^' , +a-ycos(t^) 

^ ^ kstat , kstat 

> = = H = + a — v 



stat 



_ k kstat _ fc Vfe^taf + ™2 - kstatV k^ + to' 

Vfc^ + f77,2 V A:2(^j + m2 Vk^ + m? \J k'^f^^ + m? 

e(kl,, + m2) - fc2^,,(fc2 + m2) ^ Ifc2m2 

(fc2 + m2)(fc2j^j + to2) ~(fc2+m2)2 ' 2l J- 

Note, that hi and do not depend on a and y. 



2 
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We can use this estimate in (|69|). As g" and g have uniform bounds in a and y we get 
uniform estimates for 1\ and /|: 

lull ' 



II il lU 



Jfei>-^ "2 'H "2 

Hence 

IU2 lis + 11/2 lls<Ai^'C 

with a constant C uniform in hstat- 

The integrals /f and /| can be estimated in a similar way, partial integration now be done 
with ^2 and fca 

dk,g{k) = + 1/ < , , + T for J = 1; 2 

which is uniformly bounded away from zero on the area of integration. 
So we have a uniform constant C with: 

h < fi-^C , 
and the lemma is proven for y G [a + ^,a + 1]. 

(II) Next we prove the Lemma for y < a + 1/2. 

We again have to assure, that all estimates are uniform in a and y. In the last section the 
main difficulty we had to solve was, that g' near the stationary point was increasing with k^tat 
(recall that liuiy^a+i kstat = 00). 

So on the first view it seems to be simple to have uniform estimates for y < a + 1/2 just 
by setting Q = M^. But we have to face a new problem, which is, that the stationary point 
may be very close to zero. This is problematical in the differentiation of k appearing in our 
estimates. 

For a = this problem does not appear and the lemma is also proven for y < i with a = 0. 
As the divergence only appears for small kgtat we can set kgtat < ^ (For kgtat ^ ^ the 
estimates can be done very closely to the ones of (I), setting Q = R'^). 

We solve the problem by first "cutting out" the stationary point. We split our integral: 

/= / e-'^<^^^\(\^)d?k+ [ e-'^'a'^^^xO^)d'^k^:Ii+h. 

As kgtat < 51 the stationary point is inside the ball. 

We estimate /i, writing it in spherical coordinates "centered" around the stationary point, 
by one partial integration: 



■/fci< "'1 '"i '"i 

^ -2 /" II 0X2 , ^X29^ , 0X2.9 II ,3, 

r -V/" -V/ -V/' r," 

< l|3f|+3^+3^|U.d3fc 

Jki>2k,tat '^2 "-2 "-2 

r z" z ^"2 r," 

< „-2 / 113^+3^ + 3^11 d3fc^ 
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'1 s 



< 



-B(o,\/fc;^^) 



JBio,,/!^) 9 9 9 



As X G ^ all these terms are bounded, we have 



II h \\s< Mfi-^^/k^t 

We now estimate l2- 

The first idea is to estimate this integral by two partial integrations. But the integrand still 
comes "very close" to the stationary point, where (5')^^ is not bounded. So this procedure 
will not yield uniform bound in a and y. 



The trick to get uniform bound is to redo the split (59), (60) of (I) into the integral for 
a+5 < y < a + 1. 

I2 = ll+ll+dsl{s) \s=^ 
now using fc = as the center for our Taylor-expansion. So we get: 



J]R''\B(o,\/fc;i;7) 



where 



JlR'\B(o,Vfc;;;7) 

dsl{s)\s=^ = I Xk^f{k)e^i{k)x{k)k^dkdn 



9'iO) = + a ~ y cos{-d) \k=o^ a - y cos{d) 

V fc + 

<?"(0) = dlg= , k^o^- 



Vfc' + Tn 



■3 



2 4_ 7T^2 TTL 
1 



/(k) = (5(k)~g(0)-5'(0)fc- 25"(0)fc2)fc 
^^(k) = giO)+g'iO)k + ^g"{0)k'+Cm . (70) 

As by similar argument concerning (|6^) has only one stationary point kstat ■ One can easily 
see, that 

g'(0) + g"{0)k ^a~y cos(^9) + - > a - y cos(i?) + , / „ = ^'(k) • 
Furthermore we have, that: 

g^(k) = (l-e)(g'(0)+5"(0)fc)+e5'(k). 
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It follows, that 



Therefore at k = 'k-stat (where by definition g'^{\istat) = 0) the g' has to be negative. It follows 
(recalling, that g' increases monotonously on the fci-axis), that 



< kstat < kstat ■ 

For the same reasons we have the zero point \<istat of g'(0) + g"{0)k (i.e. kgtat — ~ g"(o) )- 

< kstat < kstat ■ 

As the second derivative of g'^{kstat) is not equal to zero, we can define a function Ij^ with: 

0<M <g^:^\k -kstat \-' g^- (71) 

The integral J| includes the leading term. It can be estimated like (|6^). The other terms 
can be estimated again by partial integrations. For that we define: 



Ci := (x(k) - x(0))fc2 =: Ci^^ C2 f{k)x{k)k' =: (2^^ 

where C1.2 are bounded C°°-functions. 

We now make two partial integrations in if and three partial integrations in dsl{s) to get 
the estimates 



'1 lis 



< 



dsKs) 



< 



5, 



Ci5"(0) 



IR^s(o.v/fe;7:7) ^n3'(0)+5"(0)fc)2 (g'(0) +.9"(0)fc)3 
5fc(^9fc(^5fe%))dfcdr! lU 

5.(la.(^-M))dfcdf7|l^ 

]R^\s(o,yfc;;^) 5^ .9^ 5^ 

,„3 , 773 -3-74 173- + 3— 75-)*'='^^ 

/IR^B(o,^/fc;;:7) 5^ 5^ g^ g^ 

So we can define functions /j, j ~ 1; 5 which are bounded, with: 



)dkdVL 



2 11 < 



dkifiql + f2qf)dkdn 



9s/(s)U=dU < ^"'ll/, dk{h<ll + Uqt + hql)dkd^ 

JW\B(i).vk::7t) 



where 



91 



k-k. 



92 := 



k-k. 



28 



So it is only left to show, that duqi and dkq2 are bounded on \ J5(0, \fk^ta't)- But this is 
easy: 



dkqi = dk 



1-2 



kstat COs(l?) _|_ fc^t„t 

k 



kstat kstat C0S(1?) 



1 - 2 



fc3 



/c2 



for k > y/ kstat this term has obviously uniform bound. 

The derivative of q2 can be estimated in the same way. We only have to replace kstat by 

kstat ■ 



(III) For y > a + 1 we have no stationary point any more. So two partial integrations are 
possible without any problem. We again choose ki parallel to y 



12 lis 



< 



I '^^m'^Tw' ^''i^' '4 -"f 



x" ^x'g" x'g" xg'" 

gl2 gl3 gl3 gl3 



3^ lU d'k 

g 



(/' means dkj)- 

This integral still depends on kstat- To get an estimate uniform in kstat we use: 



5'(k) 1=2/- 



ki 



\^k^ + m? 



— a— > 1 — 
k 



= : hik) . 



It follows: 



h \\s< M"' 



x'\^x'g'\^xg"\xg"' 



(fk =: /i-^C 



4.2 Proof of equation ([T9|) 

For each k we have two eigenstates for electrons. These two eigenstates span the two dimen- 
sional spinor subspace for electrons. In the standard representation of the Dirac matrices these 
two spinor^ are: 






















fcl 




k- 











where 



fc± = fc2 ± iks Ek = Ek + m = Vk^ 



If we now take any linear combination of these spinors Sk — akS^(k) + 6(k)s^ and compute 
for example (s^,aiSk), we get (see(||)): 

■^The spinors here are not normalized! 
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(4,aisk) = ((a*(k)4*+r(k)4*),ai(a(k)4 + 6(k)4)) 

fci fe+ 

= (a2(k) + 62(k))2^feA;i . 
With the normahzation factor 



fk-\ 

ki 


\-eJ 



(4,Sk) = {a^{k) + b^{k)){El + k^) = {a^{k) + b\k}){El + 2Ekm + m^ + k^) 
= (a2(k) + b^k)){2EkiE, + m)) = (a^ik) + b^{k)){2EkE,) 



we get: 



Analogously we get: 



ki 



(•Sk: Sk) 



By linearity (|9|) follows. 



4.3 Proof of Lemma ^.4 

(a) 



To begin with, we consider the integral 



(72) 



for bounded, integrable || / \\s and j=l;2. 

For j=l it has been proven by Ikebe that I is Holder continuous. We extend this to 
i~2. Therefore we need to estimate: 

/(x + h) - /(x - h) 

for arbitrary h with h < j (We do not need to focus on ft, > i, as /(x) is bounded). We split 
the integral into: 

/(x + h) -/(x-h) = 
(t 



(73) 



B(x,v^)'l x + h-x' |2 |x 
I B{^,i)\B{^..Vh) I X + h - x' |2 I x - h - x' |2 ^^^^ ^'^ ^ 



+ 



R3\S(x,l)^ I X + 



^1 12 



X n X 



— )f{i^)d:'x'=:h+h+h. (74) 
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For 1 1 we have: 

||/i||,<2sup{||/(x)||,} / ^—^d 



So we can find a constant M < oo, so that 

II /i(x,h) ||,< mVT^ VheM^. (75) 
For I2 we have, usmg | y/h — h \ < ^Vh: 

II ^2 lis - 11/ ( | I , \ |, )/(x~x')dV||, 

JBi^.i)\B{o,Vh) I x' + h |2 I x' - h P 



x ^-h p - I x^ + h 12 1 
s(x,i)\s(o,v^) |x' + h|2|x'-h 



< sup {II /(x) iij / _ :V|2i:r ' d^^' 

Ahx' 



< sup{||/(x||,)} / I ,^.17 T^'^'^' 

xeR^ Ji3(x,i)\i3(o,yft:) I X + h rl X- h |2 

< sup{||/(x)||,} / ^dV. 

xeR'' J B{x,l)\B{Q,Vh) ^ 

So we can find a constant M < 00, so that 

II /2(x,h) ||,< mV7^ VheM^. (76) 
For I3 we have, using similar reasoning as above: 



l|/3||.< / ^||/(x-x')|UdV<8/. / ||/(x-x')||sd' 

Since / is absolutely integrable, we can find a constant M < cx), so that 

II /3(x,h) ||,< Af/i VheM^^ (77) 

We use this estimate on (|3^), observing, that multiplied by essentially of the 

form of the integrals in ([72|).Therforc: 



II ^^(x + h)-^^(x) ||,<M\/^ VheM^^ (78) 
Now we want to focus on integrals of the form ( |7^ ) for i~2 where /(x) satisfies: 

II /(x + h)-/(x) ||,<A^^/7^. (79) 



We do a similar splitting as in ( [731 ) . Now we have for Ii , using ( 79 ) : 

lUilU < 11/ I — ^(/(x' + h)-/(x'-h))dV 

JB(-^,Vh) I x-x' 1^ 

< I / I ^-—^M^/hd?x'\ . 

JB(^,Vh) I x-x' |2 
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Thus with an appropriate M < oo: 

II /i(x,h) II, < M/i VheM^^ (80) 

For we have: 



B{0,l)\B(0,Vh) 

2 



r |x'^h|^-|x' + hH 



'S(04)\B(0,\/X) 

Since the fraction under this integral is point-symmetric to zero, we can estimate the integral 
by: 



liWs < 11/ I 7 ^ h|2| T-^h 2 l(/(x~x')^/(x + xO)dV||. 

J B{o.i)\B(o,'/h) I x' - h 1^1 x' + h 1^ 

< II / I '-^7^;-'-'t"'^ I MV2?d^.' 

J B{OS}\B{0.,Vh) 



< \\ f 4 I ^ I MV2?d^x' \<\ 167rMV2 C \ ^\ d"-' 



So we can find a Af < cx) with: 

II /|(x, h) |U< M/i VheR^. (81) 
For 1-i we do the same estimations as before. 



Applying this to (35) we obtain the Holder continuity of degree 1 for c/?^. 



(b) ^ ^ 

Assume that <i5^(x) satisfies ( |3q ) and is Holder continuous of degree 1. Inserting (/3^(x) in 
the right hand side of (p^) leads to: 



F^^(x) = (i7o + 4(x))(^^(x)-y'4(x')G+(x-x')(?^(x')rfV) 

= (£;, + - (Ho + 4(x)) I 4(x')G+(x - x')^^(x')d^x' . 

For ( ^ this term has to be equal to £'fci^^(x). So we have to prove, that 

{Ho - £;fc) J G+(x - x')4(x')^^(x')rf'2;' = 4(x)^^(x) . 
In other words we have to prove, that with / Holder continuous of degree 1: 

{Ho - Ek) j G+(x - x')/(x')d3a;' = /(x') . (82) 
can be written as Q : 
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G+(x) = {Ho + Ekf— [Ho + E,)G, 
with 

(Ho - Ek){Ho + Ek)G^'' = (A - k')G^'' = S . (83) 
So for ( ^2|) we need to show, that: 

iHo-Ek){Ho+Ek) J Gf «(x-x')/(x')d3x' = (A-fc2) J Gf «(x-x')/(x')d^.T' - /(x) . (84) 
We define for e > the following function G|: 

G|.(x) G^^(x) for a; > £ G|(x) = G^^^(x)(l - ) for a; < e . (85) 

We denote 

GUx)=VGf'^ = ^fV + ^- (86) 
We split the right hand side of (^4|) into: 

(A - k') I (Gf «(x- x') - G|(x))/(x')d3x' + (A - A:^) 1 GU^) f x' . (87) 

By definition of Gj^'^ ( |8^ ) we have outside the Ball _B(0,e): 

(A-fc2)G|(x) = (A-ft2)cf«(x)=0 
So for the first summand we have: 

lim II (A - fc2) I (^Gj^G^^ _ _ x'))/(x')d=^a:' lU 



ikx 



KG 



= lim 

e^O 


II A 


/ 

/S(x 


(Gr(x-x')- 


-G|(x-x'))/(x')dV lU 


~ lim 


II V 


/ 

'B(x. 


V(Gf«(x-x' 


')-G|(x-x'))/(x')dV||. 


= lim 


II V 


/ 

'B(x. 


(V,,(Gf«(x- 

.£) 


-x')-G|(x-x')))/(x')dV lU 


= lim 


II V 


/ 

'S(x^ 


(V.-(Gr(x') 


-G|(x')))/(x-x')dV||. 


< lim 


/ ^ 

/B(x, 


II 


(V.,(Gf«(x')- 


G|(x')))^^"""'^"f" + ^""'^ 



As f is Holder continuous, the last term can be estimated by: 
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lim II (A - k') I (Gf «(x - x') - G|(x - x'))/(x')d^x' || 



< lim / I V.,(Gr(x')-GI,(x'))Af |dV 

< lim / I fcUx') - G',(x')(l - e^) - GU^')—^) M \ d'x' = 



'B{0,e) 

For the second summand, we use (pq) and the mean value theorem 



lim(A-fc2) / G|(x-x')/(x')d^a;' 



lim(A-fc2) / G|(x - x')/(x')d^a;' 



S(x,e) 



lim / (A - fc2)G|(x - x')/(x')d^a;' 



B(x,e) 



lim / (A-fc2)(e-'''=l''"'''lG|(x-x')e'*=l''-'''l)/(x')d3a;' 

-B(x,e) 



lim / (A-fc2)(g-«fe|x-x'|(;j|(^_^/)g»fe|x-x'|^j(^/-)^3^/ 
-B(x,e) 



= lim / (A + 2iA:V)(e-*'^l^-^'lG|(x-x'))e*'^'l^-^'l/(x')rf^a;' 
= lim / A(e-*'=l^-^'lG|(x-x'))e''''l^-^'l/(x')d^a;' 
= lime*'=l^-^-l/(xe) / A(e-*'^'l^-^'lG|(x-x'))dV 



'B{x,e 

where x^ G i?(x, e) using the positivity of 



Hence with Gauss' theorem and i\8 



lim(A-/c2) [ G|(x-x')/(x')rf^x' 

= /(x)lim / A(e-''=l''-'''lG|(x-x'))dV 

= /(x)lim / V(e-*'^l^-^'lG|(x-x')) -ndJ^ 

= /(x)lim / V(e-''=l^-^'lGp(x-x')) -n I x-x' |'2 rff7 

= /(x)lim/ ^i^i^.n|x-x'pdf} = /(x) 
=^0J9i3(x,e) 47r I x-x' |J 
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and (b) is proved. 

We show now, that for any k e M"^ there exists a unique solution <^^(x) of (p5[). 
Using the definition of the Ck(^) (se^ ^8|) in ( p5| ) yields: 

C£(x) = ^;k(x) - J 4(x')Gj (x - x')C^(x')d^x' (89) 

where 

t;k(x) :=- 1 4(x')Gj (x - x')¥'^(x')d3:E' . (90) 

It suffices to prove, that (^) has a unique solution for any k e M"^. For the Schrodinger 
Greens- function, this has already been proven by Ikebe We want to proceed in the same 
way. 

Let B be the Banach space of all continuous functions tending uniformly to zero as x — > cxd. 
Due to (|4^) w(x) G B. Ikebe uses the Riesz-Schauder theory of completely continuous operators 
in a Banach space [|l6j: 

If T is a completely continuous operator in B, then for any given g ^ B the equation 

f = g + Tf (91) 

has a unique solution in K if / = Tf implies that f — 0. 
Defining the integral operator T by: 

r/(x) := - 1 4(x')G+(x - x')/(x')d^x' 
and using v for g, (|9l]) is equivalent to (p9). Note, that this operator is completely continuous 



by the proof of Lemma 3.4(a) following a similar argumentation as in ]lq| Lemma 4.2. So it is 



left to show, that the integral equation 

/(x) = -J 4(x')G+(x - x')7(x')d3x' (92) 
has the unique solution / = 0. 



Obviously / = is a solution of (|92|). By virtue of (47) any solution of ( p2| ) has to be of 
order x^^. Furthermore / satisfies 

(-A-fc2+ 4)7=0 (93) 

which can be shown by direct calculation. 

Following Ikebe, / = is the only solution of 



(c) 

(c)i) follows directly from (p7|). For (c)ii) we need to work more. We exemplarily prove (c)ii) 
for j=l,2. 

Heuristically deriving (^9|) with respect to k will yield dkC- We denote the function we get 
by this formal method by . Then 

C(x) = 9fc«k(x) - J 4(x')afeG+(x - ^')a{^')d'x' - J 4(x')G+(x - ^')Cl:{^')d^x' . (94) 
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We will now show, that this integral equation has a unique solution. We define 

p(x) := 9fe«k(x)+y4(x')afcG+(x-x')Ck(x')c?V (95) 
Ck(x) C(x)-Mx) (96) 

so satisfies: 

Ck(x) = - 1 4(x')G+(x - x')p(x')d3x' - I 4(x')G+(x - x')Ck (x')d':^^' 



Since 



this integral equation again has a unique solution, so does (|9 
We will now show, that — dkC- 
We define the integral of 



i;'(x) :=- J 4(x')G+(x - x')w'(x')d^a;' G B 



k 

CksM Co(x) + / C^>sJ^)dkf . (97) 



Jo 

Obviously ^fcCk = Ck and Co = Co- Using (||) and (|5|) in (|9^) leads to 

C^(x) - Co(x)+ f Q,^(^)dk' 







= z;o(x)-y 4(x')G^(x-x')Co(x')d'a;' + y^ 9fc.«k'(x)dfc' 

- y y 4(x')afc'G+ (x - x')C^' Wx' + y 4(x')G+ (x - xOCk- {■^)d?x'dk' 
= t;k(x) - y 4(x')Gj (x - x')Co(x')d'^' 

- y y 4(x')9fc'G+ (x - x')Ck' (x')rf'^' + y ^(xOG^, (x - xOCk- {^)d?x'dk' 

- J J 4(x')9fc'(G+(x-x')Ck'(x'))'^'^'rffc' 
= «k(x) - y 4(x')G+(x - x')Ck(x')d'^' • 

So Ck satisfies (|8^). As the solution is unique, it follows, that Ck — Ck' hence 

5.Ck = Ck'- 

By (|47| ) Ck and p(x) have uniform bound, so 

sup I! 3fcCk(x) |U< oo . 

xeR3 

For the second derivative we have: 
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■ x + l x + l x + l 

The proof of the existence and uniqueness of is the same as for Sfe^", furthermore dkC is 
bounded uniformly in x. 
For dh^TT we have: 

II \\s = II ^9,(9fei;k(x)- J 4(x')5fcG+(x-x')Ck(x')rf'^') II. 

= II ^(9^k(x) - 1 4(x')5,'G+(x - x')C^(x')d^x' 



J 4(x')afeG+(x-x')9.Ck(x')rf'^') II 



Note, that 



I dlG+i^) H x'S+{^) + xdkS+{^) + dlS+{^) \< M{xk + A) . (99) 
Observing ( |99| ) and (^4|) and dkCu bounded uniformly in x, we have also, that 

-l^(a^k(x) - 1 4(x')5.G+(x - x')a..c^(x')d3x') 

is uniformly bounded in x. 

For the other summand we have get: 

sup II 1 4(x')5,2G+(x-x')a(x')d3x' lU 

< sup 11 ^ / 4(x')(x-x')a(x')d3a;' lU 

< sup II / 4(x')-4t^-;xtt(^' + i)Ck(x')rf'^' II. 

xeE3 7 (a; + l)(a;' + 1) 

< sup II / 4(x')A/(x' + l)Ck(x')d^a;' ||,< 00 



This proves (c)ii). 
(c)iii) 

The proof of (c)iii) is very similar to the proof of (c)ii). The only difference is, that we get 
new functions p(x). 

p(x) = k\^\-^Dlv^{^)) + j 4(x')fcl^l-'D^G+(x-x')Ck(x')rf'^' ■ 

To have p(x) in B one only has to assure, that k^^^~'^D'^k is bounded for | 7 |< 2, which 
follows by direct calculation. 
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(d) 

For potentials satisfying Condition A (||) the scattering system {H,Hq) is asymptotically 
complete (see Q), i.e. for any scattering state ip there exists a free outgoing asymptotic V'out 
with: 

lim II V(x,t)-Vout(x,i) IhO. (100) 

t — >oo 

We write this, using the Fourier transform V-'out of V'out: 

2 

lim ||V(x,i)-V {2n)-i^:outA^WlMd'k\\=0 . 

We shall show that 

lim II /(2^)-iv^„„4,,(k)(^£(x,i)-^£(x,i))d3fc||=0. (101) 

With that: 

lim II ^/.(x,i)-V / (2^)-'V^o„t,,(k)(^^(x,i)d3fc II 

s=l 

2 

= lim II e-^^*(V'(x) - V / i2n)-^out,s(MWki^)d^k) \\ 

s=l 

= II V(x) - / (2^)-^V^o«M(k)^kWd'fc 11= . 

s=l 

which establishes (^. For (llOlD we consider 



(2^)-5^o„t,,(k) (^^(x,t) - (^^(x, t))d^k 

(2^)-ie^^''-*V^o„t,,(k)C£(x)d3fc 

(2^)-ie^^''-*V!^o„t,3(k)<(x)d3A: 

^ j{2nrie'^''%utA^) j 4(x'^G+(x - x')Ck(x')rfVd3fc 
6(x)+6(x). (102) 



For the k- integration of we introduce (90) and (pj) and then use Lemma |3.l|, setting 



^ = a = I x - x' I; y ^ i^^x'; k' = fc; x(k') = (27r)-5^o„t^,(k') 
Furthermore we recall, that: 

a - y = 



t 
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For ^2 we set: 

M - t; a = t-i I X - x' I; y = 0; fc' - k; x(k') = (2^)-ic^(k')^out.s(k') . 
Hence by ( [l^ ) we obtain for ( |102| ) that there exists M < oo uniform in y and a, such that: 

II ei(x) +6(x) lU < Mt-i I 1 4(x')G+(x-x')(l + x')dV I 

= : Mi^iG(x) . (103) 
The integral G(x) is bounded and goes to zero in the limit a; ^ oo (see|^. This we shall use 



in the following estimate. For (101) we need to control 



lim |Ui+6lh lim ( / |Ui + 6 II' d'^;) ^ 



We split this integral into three parts, which are time dependent by introducing for all e > 0: 
^^(x) = ^B{o,et) 1 the indicator function of the set _B(0, ei) 

Pe(x) = lB(0,t)\S(0,et) 
Pout(x) = lR3\s(0,t) 

thus splitting our integral into: 

hm / II a + 6 \\l d?x = hm / p,(x) || + 6 ||^ d?x 

+ lim / p,(x) II 6+6 lis <i^x 

+ lim / pout(x) II 6+6 lis d^x 

= : h+h+h. (104) 

The last part of this integral is the part, that lies outside the light cone. For large times, all 
wavefunctions which are solutions of the free or full Dirac equation will lie inside the light cone. 
By virtue of (pi): 



by ([0^: 



hm II pout(x)'0out,s(x) II = 
lim II pout(x) / {27r)-^outA'^)'pU^)d^k \\ = 



lim II pout(x) / i2Tr)- HoutA^)vU^)d''k |h 



By ( |102| ) it follows, that: 



/3 = lim II pout(x) / ^out.sO^) (^^(x) - ^^(x)) d^k 11= 
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Now we use (1103) on 



h < lim (sup{G(x)})2f-3 — (et)3 = Ce^ 

^ X<£t 'J 



Since e is arbitrary, = 0. 
For /2 we have: 



I2 = lim 



lim I / M^)U,+^2\\W^\ 

t^oo J 

lim I Af2 / t-3^e(-x)G'(x)2d3^ I 

:^00 J 



t 

< lim sup I G(x)2 1= 

t— 00 ^>gj 

and (^) is proved. 

We first prove (^0|) for wavefunctions, where ipout is in D L'^. The general result can then 
be obtained by density arguments. 

Therefore we again use the unitarity of the time propagator: 



(27r)-S(^«(x)>(x))d^a; = lim / (27r)-Se^^* (^^(x), e"^^ V(x))rf 



t~*oo 



X 



hm e'^* /(2^)-2(^^(x),V;(x,t))d^a; 

lim e'-^* / (27r)-i(^^(x),V(x,i))d3x 
*^°° Jb(o,r.) 

+ lim e^^* / (2^)-i((^^(x),i^(x,t))d3a; . 



By asymptotical completeness (|100|) we obtain therefore 



(2^)-5(^^(x),V(x))d3x = hme*^*/ (27r)-5 (^^(x), V',„t(x, i))^^^ 



B{0,R) 

Et / /~s /„\ /„ j-N\j3„ 



+ hm e'^* / (2^)-^(<?^(x),^out(x,i))d^2: . 

JR3\B(0,_R.) 

By the free scattering into cones theorem, the first integral of the right hand side goes to zero 
because any freely evolving wavefunction leaves any bounded region. For the second integral 
we write for all i? > 0: 



(27r)-5(^^(x),^(x))d3:E = lim lim e'^* / (27r)-5 ((^^(x), Vw(x, t))d^:r 



i?,^oo t^oo 



3\S(0,i?) 



lim lime^^*/ (27r)-i (C^(x), Vout(x, i))^^^ . 



3\B(0,i?) 



Using Lemma 3.4 (c)i), the second integral on the right hand side becomes: 
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lim lime*^*/ (27r)-3 (C^(x), V-outlx, i))^^^ | 



< lim ^11 lime*^* / (27r)-i V^out(x, ^d^x || 

R^oo K t^oo JwL^\B{Q,R) 

= lim ^11/ (2^)-iV^out(x, 0)^3:^11=0. 



Therefore: 



(27r)"*(^^(x),V'(x))d3a, 



= limlime'^M (27r)- ^ (<^^(x)>„,t(x, 

fl^oo t^oo Jm.^\B(0,R) 

= lim lim e^^*/ (27r)-i ((^^(x), Vout(x, i))^^^ 



R3\B(0,i?,) 



= lim lim e'^M (27r)-3(<^^(x),^out(x,i))d'^x 

R-^00 t^oo 

= lim / (27r)-5(^^(x),V;out(x,0))d3a; = Vo«t,5(k) . 



and ^dj ) is proved. 



4.4 Proof of Lemma 2.2 



First we want to prove 

Let V^outCk) G tj. According to we have for any n £ Nq: 



2 . 

iI"V(x) = ^ /(2^)-iff"<^^(x)V^„„t,,(k)d3fc 

s=l 

Since V'out(k) decays faster than any polynom, this term is bounded and in i^(^C"' for all 
n G No. As the potential 4 € C°° , also 

{H~4- /3m)"V(x) = TV'(x) 

is bounded and in C"' for all n G Nq. 
Furthermore we have, using (pq) in (B9|): 
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s=l 
2 

= ^ / (27r)-i^^(x)i?,"^o„t,.(k)d3fc 

2 

/ (2^)-i f 4{^')G+{^-^')(pi{^')d^x'E^iPoutA^)d'k =: /i + /a . 

s=l 

/i is the Fourier transform of EJ!'tpout,s{'k). As i?^V'o«t,s(k) € G, h lies in G. 
Next we write for J2: 



where 



^^(k,x,x') e'^^(|x-x'|-^) '^fc(x • (105) 



X — X' 

We make now two partial integrations under the k-integral, which is possible by Fubinis theo- 
rem: 



h = i'^^y'^ J 4{^')e:''^F{k,^,^')d''x'E^k^i^out.A^)dkdn 

= -y"— /(27r)-M / 4(x')e'""9,2(^^(k,x,x')dV-B^fc2V^o„t,,(k))dMf7 
s=i ^ ■> -'0 J 

= -E^ / ^°°/(27r)-i4(x')e^'"a^^(k,x,x')£;,"fc^V^o«t,.(k) 

+ 2dkF{k, X, x')5fc {E]:k^iPout,s (k)) dfcdf^d^x' 

/(2^)"^ / 4(x')e''"i^(k,x,x')dVa^(i;fc"A;2v^o.t,.(k))dM17 



For /4 we can write, using the definition of F (105) and (p. 



2 . -, 
x2/4 = ^ /(27r)-i^^(x)a2(i?;'A:Vw.(k))yjd^fc 

2 

i2nrhU^)dUEj:k%^tA^))^d'k 
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As V'out e G, dKEJ^P^outA^)) w l^^s ^^'^ '^^^^ x^d]^h for n € No- 
Under the k-integral in /3 one more partial integration is possible. 

2 

/3 = - E 4 / (2^^)"^ / 4(x')i^(k,x,x')d3x 

where 

i?(k, X, x') := 5fc {dlF{]^, X, ^)Eie^ouUs (k) + 2afc^^(k, x, x')afc (£;^fc2^o„t,,(k)) ) . 
Due to Lemma (c) S^c^klx'-* < Afa;'. Furthermore we have, that 

I ^fce^'^'^l^-^'l-^) 1 = 1 (I x-x' I -a;)e*'^(l''-^'l-^) |< .t' | e^'^d^-^'l"^) | . 
It follows, that (remember the definition of F ( |105D ) 

II i?(x,x) ||.<M2 , , , . 



So due to (^7|), with Condition B (^) on the potential, the integral 

j 4(x')f^(k,x,x')d^a;' 

decays as fast as or faster than x^^, so x^Iz is bounded. It follows, that x^Iz hes in for 
n e No. 

The proof, that xd^ip € is similar as above, just with one partial integration less. It 
follows, that tp G G- 
It is left to prove 
By Lemma |3.4| (b) it follows, that 



(2^)-i(^^(x),(i/o + 4)V'(x))d3x. 



For ip & G, the right hand side is integrable, so i?fe-0out.s(k) is bounded. As 4 G , this can 
be repeated, so i?^i/'o«t,s(k) is bounded for any n e N. 
Since — \/k^ + m^ > fc, it follows, that 



A:"V'o«t,s(k) < oo . 



Equivalently we get: 



EJ^k^^^-^Dl^outA^) = [{2TT)-i{k\^\-^Dl^l{K),H"^P{x))d^x 



With (c) of Lemma 3.4 it follows, that for tp G G these terms are bounded for j — 1,2, n £ No 



and I 7 l< 2. So ipout^sO^) G G- 
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